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!n  this  article  we  consider  two  versions  of  two-on-two  homogeneous  stochastic 
combat  and  develop  expressions,  in  each  case,  for  the  state  probabilities.  The 
models  are  natural  generalizations  of  the  exponential  Lanchester  square  law 
model.  In  the  first  version,  a  marksman  whose  target  is  killed  resumes  afresh  the 
killing  process  on  a  surviving  target;  in  the  second  version,  the  marksman  whose 
target  is  killed  merely  uses  up  his  remaining  time  to  a  kill  on  a  surviving  target. 
Using  the  state  probabilities  we  then  compute  such  important  combat  measures 
as  (1)  the  mean  and  variance  of  the  number  of  survivors  as  they  vary  with  time 
for  each  of  the  sides.  (2)  the  win  probabilities  for  each  of  the  sides,  and  (3)  the 
mean  and  variance  of  the  battle  duration  time.  As  an  application,  computations 
were  made  for  the  specific  case  of  a  gamma  (2)  interfiring  time  random  variable 
for  each  side  and  the  above  combat  measures  were  compared  with  the  ap¬ 
propriate  exponential  and  deterministic  Lanchester  square  law  approximations. 
The  latter  two  are  shown  to  be  very  poor  approximations  in  this  case. 


1.  INTRODUCTION 

The  principal  motivation  for  this  work  is  the  development  of  more  realistic 
small-to-moderate-size  firefight  models.  It  is  an  extension  of  the  work  started 
in  references  [1]  and  [7]  which  treat  the  one-on-one  and  homogeneous 
two-on-one  stochastic  combat  models,  respectively. 

The  overall  framework  within  which  all  these  works  lie  is  described  in 
reference  [3]  where  (1)  a  comprehensive  examination  is  presented  of  the 
nature  of  combat  and  the  status  of  corresponding  theory,  and  (2)  a  proposal 
toward  a  theory  of  combat  is  set  forth.  The  basic  conclusion  reached  in 
reference  [3]  is  that  current  modeling,  no  matter  how  “realistic”  it  is  claimed 
to  be,  is  deficient  and  not  based  on  any  firmly  established  theory.  One  only 
needs  to  look  at  an  application  of  the  classical  Lanchester  square  law  to  large 
numbers  of  opposing  forces  (inserted  as  initial  numbers)  to  see  the  following 
two  egregious  assumptions,  concerning  the  nature  of  combat: 


*This  work  was  partly  supported  by  the  Army  Research  Office,  Contract  No. 
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/.  Coinhcil  is  homogeneous.  This  in  effect  assumes  that  in  lOOO-on-l ()()() 
battles  all  the  forces  are  simultaneously  engaging  the  enemy  in  the  same  way. 
Common  sense  dictates  that  the  effects  of  cover,  concealment,  weather, 
terrain  obstacles,  terrain  corridors,  the  effectise  range  of  different  weapons, 
countless  other  environmental  factors,  and  tactical  deployment  force  the 
opponents  into  many  smaller  firefights.  An  excellent  field  study  in  this  con¬ 
nection  appears  in  reference  [9]  where  it  was  found  that  a  structured  relation¬ 
ship  for  the  minibattles  may  be  constructed  showing  how  some  battles  occur 
simultaneously  (or  in  parallel)  while  others  involve  the  participants  in  a  series 
of  battles  through  time. 

2.  Combat  is  deterministic.  To  believe  that  combat  is  not  stochastic  would 
require  that  if  a  firefight  could  be  repeated  under  exactly  the  same  conditions, 
every  event  that  occurred  in  the  first  replication  would  reoccur  on,  say,  the 
second  replication  in  exactly  the  same  order  at  exactly  the  same  time.  This 
truly  would  be  incredible.  Stochasticity  in  combat  is  indeed  significant  and  an 
extensive  discussion  of  it  may  be  found  in  reference  [2].  In  addition,  reference 
[9]  describes  an  excellent  experiment  showing  not  only  stochasticity  in  a 
combat  process  but  also  its  considerable  variance. 

We  hold  to  the  view  described  in  reference  [3]  that  the  ultimate  develop¬ 
ment  of  realistic  combat  models  involving  large  numbers  of  weapons  will 
depend  on  successfully  modeling  (1)  the  decomposition  of  the  large  battle  into 
the  separate  small  engagements  (in  this  connection  see  references  [4]  and  [5]), 
and  (2)  the  attrition  process  in  the  separate  engagements.  This  study  is 
directed  toward  (2),  the  realistic  modeling  of  small-to-moderate-size  engage¬ 
ments. 

Furthermore,  we  believe  that  to  achieve  realism  in  the  attrition  process  we 
must  proceed,  as  it  has  generally  been  done  in  the  physical  sciences,  from  the 
simple  to  the  more  complex,  and  so  this  study  is  the  obvious  successor  model 
to  those  described  in  references  [1]  and  [7],  It  provides  what  we  believe  to  be 
substantially  more  realism  in  the  extant  two-on-two  stochastic  (more  ac¬ 
curately,  exponential)  Lanchester  square  law  combat  models.  It  does  so  by 
removing  the  extremely  simplifying  assumption  that  the  interkilling  time 
random  variable  is  the  same  negative  exponential  distribution  (NED)  from  kill 
to  kill  (this  is  a  consequence  of  assuming  that  both  the  single-shot  kill 
probability  and  the  NED  interfiring  time  random  variable  are  the  same  from 
round  to  round).  Allowing  both  the  single-shot  kill  probability  and  interfiring 
time  random  variable  to  vary  from  round  to  round,  of  course,  complicates  the 
analysis  substantially  and  numerical  techniques  are  required  to  produce  state 
probabilities  in  any  specific  case.  We  do,  however,  retain  the  homogeneity 
assumption  that  all  combatants  on  a  side  possess  identical  characteristics  (but 
which  are  not  necessarily  the  same  for  both  sides). 

The  two-on-two  raises  the  question  of  how  a  marksman  whose  original 
target  is  killed  by  another  handles  a  surviving  target.  In  this  article  we 
consider  two  versions  of  how  this  situation  is  dealt  with.  In  the  first  version  the 
marksman  starts  the  killing  process  all  over  again  on  the  surviving  target,  and 
in  the  second  version  the  marksman  uses  up  the  remaining  time  to  a  firing  (or 
a  killing)  on  the  surviving  target.  In  the  exponential  Lanchester  ca.se  there  is 
no  difference  in  these  two  versions  because  of  the  no-memory  property 
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associated  with  the  exponential  random  variable.  In  each  of  these  lersinns  we 
have  allowed  all  possible  values  for  a  side's  breakpoint  (the  number  of  survivors 
at  the  time  the  side  loses). 

In  the  following  sections  the  problem  is  precisely  formulated  for  both 
versions  and  the  state  probabilities  are  then  derived.  These  state  probabilities, 
which  are  functions  of  time,  contain  all  the  pertinent  information  about  this 
nonstationary  stochastic  process,  and.  in  fact,  in  terms  of  them  we  can  write 
expressions  for  the  following  important  combat  measures  of  effectiveness; 


1.  The  mean  and  variance  of  the  number  of  survivors  as  they  vary  with  time  for  each  of 
the  sides.  At  /  =  rhi.s  jives  the  mean  and  variance  of  the  number  of  survivors  of  the 
battle  for  each  of  the  sides. 

2,  The  win  probabilities  for  each  of  the  sides. 

.1  The  mean  and  variance  of  the  battle  duration  lime  To  Actually  one  may  write  the 
distribution  function  of  battle  duration  Ftd'''-  terms  of  the  state  probabilities  and 
the  equivalent  probability  density  function  /tdI''  in  lerms  of  their  derivatives. 
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The  technique  used  in  this  article  results  in  n-fold  iterated  integrals  for  the 
state  probabilities.  The  integrands  are  complex  products  of  the  density  and 
complementary  distribution  functions  of  the  interkilling  time  random  variables 
of  both  sides.  The  dimension  of  the  ii-fold  integrals,  for  any  particular  state 
probability,  is  equal  to  the  number  of  kills  corresponding  to  the  state.  Thus, 
e.g.,  in  the  zero-breakpoint-for-both-sides  case  we  get  up  to  three-dimensional 
iterated  integrals;  and  the  computation  time  to  get  all  state  probability 
functions  plus  eight  overall  battle  parameters  is  from  two  to  three  hours  on  a 
Sperry  1 100/82  when  the  interfiring  time  random  variables  are  gamma(2)  for 
both  sides. 

We  also  note  here  that  the  present  method  certainly  can  be  extended  to  the 
three-on-two  and  even  the  heterogeneous  versions  of  battles  up  to  that  size, 
where  by  heterogeneity  we  mean  that  the  characteristics  of  the  combatants  on 
a  side  arc  allowed  to  be  different;  and  these  extensions  are  presently  under 
way.  Desirable  as  it  is  to  have  exact  solutions,  in  view  of  the  computer  times 
involved  in  using  the  present  approach,  we  feel  that  alternative  methods  must 
be  considered.  These  include  the  following; 

1.  Simulation.  We  have  already  developed  arbitrary  fiu-on-bi,  versions  of 
these  models,  including  the  requisite  statistical  techniques  for  their  use.  We  are 
presently  in  the  evaluation  process. 

2.  Approximations.  We  have  defined  a  nonhomogcncous  Poisson  process 
approximation  which  results  in  vastly  simpler  analytical  expressions  ft'r  all  the 
parameters  of  interest  and  for  which  the  computer  time  required  to  produce 
answers  in  any  specific  case  is  substantially  less  than  required  by  our  exact 
model.  Furthermore,  the  simulation  version  of  this  approximation  is  much 
faster  than  the  exact  version.  We  are  presently  evaluating  this  approximation. 

3.  Other  e.xact  techniques.  These  may  result  in  simpler  expressions  for  the 
state  probability  functions  than  we  have  developed  and,  therefore,  presumably 
require  much  less  computer  time  in  any  specific  applicatiiin.  As  of  this  time  we 
have  not  come  up  with  any  sueh  technique  nor  have  we  been  successful  in 
simplifying  our  results. 
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2.  THE  MODELS 

Two  sides.  A  and  B.  conduct  a  continuous  engagement  satisfying  the 
following  essentially  Lanchester  square  law  assumptions; 


1 .  There  are  initially  two  on  each  side. 

2,  Hverv  niemher  of  A  side  picks  a  B  opponent  at  random  (all  arc  visible  and  in  range). 

.V  Each  marksman  fires  until  killed  or  until  his  target  is  killed  and  resumes  firing 

immediately  on  the  survivor  in  one  of  two  distinct  ways  which  are  described  below  and 
denoted  as  Versions  I  and  2. 

4.  Tbe  iiilerkilling  time  random  variable  does  not  change  from  kill  to  kill  and  is  identical 
for  all  members  of  the  A  side. 

5.  All  lire  independently. 

(i.  The  ammunition  supply  is  unlimited. 

7.  Simular  assumptions  apply  ti>  the  B  side. 

S.  The  battle  continues  until  one  side  reaches  its  breakpoint  (the  number  on  a  side  at  the 
time  it  loses). 

The  two  modes  of  resuming  firing  on  a  survivor  are  as  follows. 

Version  I.  Consider  a  given  marksman  firing  at  a  target.  Whether  his 
target  is  killed  by  him  or  the  other  member  of  his  side,  he  resumes  afresh  the 
intcrkilling  process  on  the  survivor. 

Version  2.  If  the  marksman's  target  is  killed  by  him,  he  starts  afresh  the 
interkilling  process  on  the  survivor.  If.  on  the  other  hand,  his  target  is  killed  by 
the  other  member  on  the  side,  his  remaining  time  to  a  firing  (or  a  killing)  is 
carried  over  to  the  survivor.  The  jargon  we  use  to  describe  Versions  1  and  2 
are  “reselect  on”  and  “reselect  off.”  respectively. 

When  we  consider  all  the  possible  breakpoints  for  sides  A  and  B.  respec¬ 
tively,  we  get  a  total  of  five  models.  These  are  shown  in  Table  1  below  with  a 
model  numbering  scheme  for  ease  of  reference. 

3.  GENERAL  SOLUTIONS 

Our  solution  technique  will  depend  on  knowing  the  interkilling  time  random 
variable’s  density  and  complementary  distribution  functions  for  each  side.  So 
long  as  they  can  be  described  analytically  or  developed  in  tabular  form  (with 
exact  entries  using  some  numerical  techniques  or  estimated  using  Monte 
Carlo)  the  formulas  derived  in  this  article  may  be  used  to  compute  the  state 


Table  1.  A  brief  description  of  the  five  rr  dels. 


Breakpoints 

Model  no. 

Side  A 

Side  B 

"Resclcct” 

1.1 

0 

0 

on 

1.2 

0 

0 

off 

2.1 

I 

0 

on 

2.2 

1 

0 

.)ff 

3 

1 

1 

Not  material 

Gafarian  and  Manion;  Stochastic  Combats 


725 


probabilities  and  subsequently  overall  combat  and  time-varying  characteristics 
(this  means,  for  example,  that  the  models  handle  variable  single-shot  kill 
probabilities  and,  for  that  matter,  variable  interfiring  time  random  variables). 
However,  questions  regarding  ammunition  consumption  cannot  be  addressed 
because  our  solution  technique  perforce  loses  all  information  on  number  of 
rounds  fired. 

In  this  article  the  solutions  to  Models  1.1  and  1.2  are  given  in  some  detail. 
These  two  have  the  largest  number  of  states  and  are  the  most  difficult  to  treat. 
In  particular,  we  begin  this  section  with  Model  I.l  and  proceed  far  enough 
along  to  give  the  reader  some  notion  of  how  state  probabilities  are  derived  and 
then  place  the  remainder  of  the  derivations  for  Model  1.1  and  all  of  Model  1.2 
in  Appendix  1.  There  we  present  also  only  the  results  for  the  other  breakpoint 
combinations  (three  more  models). 

Notation  we  use  throughout  the  remainder  of  this  article  is  as  follows: 


Ufl  =  the  initial  number  on  side  A  (at  time  0) 

Of  =  breakpoint  for  side  A.  i.e.,  the  number  on  side  A  at  the  time  the  A  side 
loses  (breaks  and  runs) 
bo  =  the  initial  number  on  side  B  (at  time  0) 

bf  ~  breakpoint  for  side  B,  i.e..  the  number  on  side  B  at  the  time  the  B  side 
loses  (breaks  and  runs) 

/a(<)'  Ga(')  =  density  function  and  complementary  distribution  function  for  the  time-to- 
kill-of-a-passive-target  random  variable,  side  A 
/aft).  Gaft)  =  density  function  and  complementary  distribution  function  for  the  time-to- 
kill-of-a-passive-target  random  variable,  side  B 
A(r)  =  random  variable,  number  alive  on  side  A  at  time  / 

B(i)  =  random  variable,  number  alive  on  side  B  at  lime  i 
Pabil)  =  P(A(/)  =  a.  B{l)  =  *].  a  stale  probability  function 
m^U)  =  E[A(i)],  expected  value  of  A(() 
mgii)  =  £(B(/)j.  expected  value  of  B(i) 

O’Alt)  =  standard  deviation  of  A(/) 

(Tgii)  =  standard  deviation  of  B(l) 

P[i]  =  probability  i  side  wins,  i  =  A.  B 
Tq  =  random  variable,  lime  duration  of  combat 
Gto^i)  -  complementary  distribution  function  for 
fir„  -  expected  value  of  Tp 
ir-To  =  standard  deviation  of  Tp 

V,  =  mean  interkilling  lime  on  side  /.  i~  A.  B 

tA  =  •/"a  =  r' s  kill  rate  (attrition  coefficient  for  side  B), 

fn  =  I/«'b  =  B's  kill  rate  (attrition  coefficient  for  side  A) 


and  whenever  the  single-shot  kill  probability  and  interfiring  time  random 
variables  are  the  same  from  round  to  round  we  use  the  notation 


Pi  =  the  constant  kill  probability  of  all  constants  on  side  i,  i  =  A,  B 
p,  =  mean  interfiring  time  on  side  i,  i  =  A,  B 
'a  =  '/"a  =  Pa/Ma 
fB~  ^l‘'o  ~  Pal  Pa¬ 


in  some  of  our  calculations  we  use  the  backward  recurrence  time  technique 
to  write  the  state  probability  equations.  If  at  time  t  we  define  Y  to  be  the  time 
since  the  last  event  (kill),  then  the  first-order  probability  that  an  A  marksman 
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will  kill  in  the  interval  {t,  t  +  A)  is  given  by 

rAy)A  =  fA(y)/GAy)^. 

and 

rB(y)A  =  /B(y)/GB(y)A 

for  a  B  marksman.  The  r{y)s  are  the  instantaneous  kill  rates  for  each 
marksman.  See  [1]  for  a  discussion  of  the  backward  recurrence  time  technique. 

.Model  1.1  (flo  =  2,  a,  =  0;  bo  =  2,  bf  =  0;  “Reselect  On”) 

The  a;ijlyses  for  both  Models  1.1  and  1.2  proceed  in  the  order  shown  in 
Figure  1 ;  namely,  we  will  write  an  expression  for  p2:(f).  then  pizd)  and  (0. 
etc.  In  many  of  our  considerations  it  becomes  convenient  to  delineate  the 
various  aiming  configurations.  Because  of  the  homogeneity  of  the  combatants 
on  each  side  the  total  number  of  distinct  configurations  is  four,  which  is  the 
product  of  the  two  ways  the  As  aim  at  the  Bs  (both  As  aim  at  the  same  B  or  at 
different  Bs)  and  the  two  ways  the  Bs  aim  at  the  As  (both  Bs  aim  at  the  same 
A  or  at  different  As).  The  analysis  for  Model  1.1  is  carried  out  in  terms  of  the 
initial  aiming  configurations.  Thus  what  we  do  is  break  up  each  of  the  states 
shown  in  Figure  1,  except  (2.  2),  into  subsets  that  are  associated  with  the  initial 
aiming  configuration.  We  define  these  subsets  now  and  the  associated  state 
probabilities. 


Figure  1.  Sequence  of  states  in  the  two-on-two  combat  with  zero  breakpoints. 
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(1.2) 

Figure  2.  Decomposition  of  state  (1,2). 


Let  us  first  consider  the  state  (1,2)  and  break  it  up  into  states  (1,2)'"  and 
( 1 . 2)''*  as  shown  in  Figure  2.  State  (1 , 2)'"  is  the  state  (1,2)  achieved  from  the 
initial  aiming  configuration  in  which  both  Bs  are  aiming  at  the  same  A  and 
( 1 .  2)'''  is  the  state  (1,2)  achieved  from  an  initial  aiming  configuration  in  which 
each  B  is  aiming  at  a  different  A.  The  corresponding 

/>i’(f)  =  pV’it) + 

Similarly, 

P:i(f)  =  P;^i’(»)  +  p-i’i’lt). 


where  the  state  (2,  !)'”  is  state  (2,  1)  achieved  from  an  initial  aiming 
configuration  in  which  both  As  are  aiming  at  the  same  B  and  state  (2,  I)'"”  is 
state  (2,  1)  achieved  from  an  initial  aiming  configuration  in  which  each  A  is 
aiming  at  a  different  B.  The  reason  for  naming  these  states  as  (2,  !)'■’’  and 
(2,  1  )'■"  instead  of  (2,  1 )'"  and  (2,  I)'"’,  respectively,  will  become  clear  when  we 
write  down  the  decomposition  for  the  states  (1,  I),  (0,  1),  and  (1,0). 

As  far  as  states  (0,2)  and  (2,0)  are  concerned,  they  may  be  decomposed 
into  the  states  shown  in  Figure  .7,  in  which  we  have  that 


,  ,  O),  ,  ,  (21,  , 

Pii’f  f)  =  Pi)2(  t)  +  P(i;(f). 


and 


P2ii(')  =  P2a(')  +  pil(t). 


As  above,  the  superscripts  (1)  and  (2)  correspond  to  both  Bs  aiming  at  the 
same  A  and  each  B  aiming  at  a  different  A,  respectively;  whereas,  the 
superscripts  (3)  and  (4)  correspond  to  both  As  aiming  at  the  same  B  and  each 
A  aiming  at  a  different  B,  respectively. 


(0,2)  (2,0) 


Figure  3.  Decomposition  of  states  (0,  2)  and  (2,  0). 
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f  {o,n‘^'  \ 

[  (1.1)'^’  \ 

1  (1,0)'^’  ^ 

V  (0,1)'^'  / 

1  (1.1)*^'  / 

V  (1.0)*^'  J 

V(0.l)“*_^ 

V|1.0)^ 

(O.t) 

(1.1) 

(1,0) 

Figure  4.  Decomposition  of  states  (0, 1),  (1. 1),  and  (1, 0). 


We  finally  decompose  the  states  (1,1),  (0,  1),  and  (1.0)  into  the  states  shown 
in  Figure  4.  In  thinking  of  these  states  the  reader  should  clearly  keep  in  mind 
the  meaning  of  the  superscripts  (/),  /  =  1, 2,  3, 4.  They  are 

(1)  Initially  both  Bs  are  aiming  at  the  same  A  and  the  first  kill  is  of  an  A  by  a  B. 

(2)  Initially  each  B  is  aiming  at  a  different  A  and  the  first  kill  is  of  an  A  by  a  B. 

(3)  Initially  both  As  are  aiming  at  the  same  B  and  the  first  kill  is  of  a  S  by  an  A. 

(4)  Initially  each  A  is  aiming  at  a  different  B  and  the  first  kill  is  of  a  B  by  an  A. 

Bearing  this  in  mind,  (0,  1)"’,  say,  is  a  state  resulting  from  a  B  achieving  the 
first  kill,  follow  d  by  the  surviving  /\  killing  a  B.  and  finally  the  surviving  B 
killing  the  surviving  A,  and  all  this  from  an  initial  aiming  configuration  in 
which  both  Bs  were  aiming  at  the  same  A. 

/.  P2:(f)-  We  begin  by  setting  down  immediately 

P22(t)  =  (C,,(t)f{Gait))-.  (1) 

which  merely  states  that  each  of  the  contestants  has  a  time  to  kill  >  f. 

2.  PniO,  P2i(0-  To  compute  p^Cf)  consider  Figure  5  below,  which  shows 
that  an  A  is  killed  by  a  B  in  the  time  interval  (t-  ■q-  drt,t-  rf)  with  no 
subsequent  killings  until  beyond  t.  Now  either  both  Bs  were  aiming  at  the 
same  A  or  they  were  not.  Each  of  these  initial  aiming  configurations  has 
probability  of  We  now  define  p' '2 ( f,  tj)  drj  =  probability  that  both  Bs  are 
aiming  at  the  same  A,  one  of  the  Bs  kills  an  A  in  (f  -  17  -  dr},  f  -  17),  and  there 
are  no  other  killings  until  beyond  t,  and  pi2(f,  17)  dq  =  probability  that  each  B 
is  aiming  at  a  different  A,  one  of  the  Bs  kills  an  A  m  {t  -  q  -  dq,  t  -  q),  and 
there  are  no  other  killings  until  beyond  i. 

Once  we  write  expressions  for  pV?!!.  v)  ^^d  puif.  tj)  we  can  get  pnit)  ^s 

Pi2(f )  =  PV2V)  +  Pn(l)  =  {  dq  {pait,  q)  +  Pu(f,  Tj)). 

Jo 


X  a  B  kills  an  A  in  (t-ri-arj,  t-T]) 

Figure  5.  Definition  of  the  variable  17  for  the  computation  of  pn(i). 
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Let  us  first  consider  pV:(t.  T)WiT7.  It  can  be  written,  because  of  the  in¬ 
dependence  assumption  of  all  the  fircrs,  as  5  times  the  product  of 

2/n(t  ~  =  pmbabilily  that  one  of  the  two  Bs  kills  an  A  in  the  time  interval 

1(  -  dr).  I  -  Tj) 

Cin(i  -  Tj)  =  probability  that  nonkilling  B  has  a  kill  time  >  1  -  r) 

=  probability  both  Bs  start  over  again  on  surviving  A  at  time  ( -  tj  and 
have  a  kill  time  >  rj 

GaI'  -  rj)  =  probability  that  killed  A  had  a  time  to  kill  >  1-  rj. 
and  finally  GaIi)  =  probability  that  the  surviving  A  has  a  time  to  kill  >  t.  Thus 
p\2( L  n)  =  ( 1  /2)[2/b( t-r))GB(t-vHGs( r,))’  f  -  t,)  Ga( f )]• 

Similarly, 

pn(t,  v)  =  (l/2)[2/B(t-v)Gs(v)GB(i}GA(l-v)GA(l)]. 

where 

2/b(<  “  drj  =  probability  one  of  the  two  Bs  kills  an  A  in  the  time  interval  (r-rj- 

dij.  I  -  T)) 

Gfifr?)  =  probability  that  the  killing  B  starts  the  firing  process  all  over  on  surviving 
A  and  has  a  time  to  kill  >  77 
Gfl(r)  =  probability  nonkilling  B  has  a  time  to  kill  >/ 
yp-  probability  that  killed  A  had  a  time  to  kill  >1-  t). 

and  finally  G^ff)  =  probability  that  the  surviving  A  has  a  time  to  kill  >  I.  Thus 
we  get 

p,2(f)  =  GaIO  [  d-nfnit-  r\)GB(l-  b)(Gb{ti))'Ga(i  ~  17) 

Jil 

+  GB(f)GA(f)[  dr}fB(i- ri)GB(ri)GA(t- T^).  (2) 

Jtt 

Obviously  by  an  interchange  of  subscripts  we  may  write 

P2i(0  =  Gn(l)  [  dT}/A(/- t))Ga((- T))(GA(T)))^G8(f- tj) 

Jo 

-I- GB(/)GA(f)  [  drtfAii  -  v)GA('n^G„(l  -  ri).  (3) 

The  remainders  of  all  state  probability  derivations  appear  in  Appendix  1. 


Combat  Figures  of  Merit 

The  transient  state  probabilities  p;2(f),  pi:(<).  P2i(0<  and  pn(i)  along  with 
the  absorbing  state  probabilities  p:„(l).  p„2(i).  Pind).  and  p„i(l}  provide  all  the 
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information  necessary  to  compute  the  following  commonly  used  figures  of 
merit. 

1 .  The  expected  value  and  standard  deviation  of  the  survivors  on  both  side 
A  and  side  B  as  a  function  of  time.  For  side  A 

nu(/)=  lP[A(r)=  l]  +  2P[A(t)  =  2] 

=  l(p!o(0  +  Pii(t)  +  Pi2(f))  +  2(p2„(r)  +  p:,(r)  +  p:2(0).  (4) 


Computing  £[A^(r)]=  lP[A(r)  =  l]  +  4P[A(t)  =  2] 
will  then  provide 


aAt)  =  iE[A\t)]-m\(t))'‘\ 


Similarly, 


msit)  =  l(pOl(t)  +  Pll(t)  +  P2l(r))  +  2(po2(t)  +  Pl2(0  +  P22(0 


and 


(5) 

(6) 


o-B(r)  =  (£[B'(t)]-m|(t))'^  (7) 

where 

£[B^(r)]  =  1  P[B(t)  =  1]  +  4P[B(t)  =  2]. 

In  particular  we  get  the  expected  values  and  standard  deviations  of  the 
number  of  survivors  on  the  A  side  and  B  side  by  letting  r^oo  in  Eqs.  (4)-(7). 

2.  The  expected  value  and  standard  deviation  of  To  the  time  duration  of 
combat.  These  are  computed  using  the  well-known  integral  formulas,  derived 
by  an  integration  by  parts,  for  the  first  and  second  moments  of  a  generic 
non-negative  random  variable  X  with  finite  second  moment,  and  com¬ 
plementary  distribution  function  G(x),  namely 

£[X]=f  Gix)dx,  £[X']  =  2f  xG{x)  dx. 

Jk  Ji) 

The  standard  deviation  is  computed  as  (£[X^]-  £^[X])'^^.  In  the  case  of  time 
duration  of  combat  we  use  the  obvious  result  that 


Thus 


and 


Gto(0  =  P22(0  +  P2l(f)  +  Pl2(0  +  Pll(0. 


=j: 


(P22(f)  +  P2l(0  +  P2l(t)  +  PmU)  +  Pi  l(t))  dt 


(8) 


(9) 
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where 


E[Th]  =  2\  tip22(l)  + P2\{t)  + P\2(t)  + pn{t))  dt.  (10) 

Jo 

3.  The  probabilities  of  win  by  the  A  side  and  B  side.  These  two  prob¬ 
abilities.  P[A]  and  P[B],  respectively,  are  obviously  given  by 

P[A]  =  lim  (p;„(f)  +  p,,)(t))  .-I) 

and 

P[B]=  1  -P[.A]  =  lim(p„.(r)  +  poi(f)).  (12) 


4.  COMPARISONS  BETWEEN  SOME  LANCHESTER  MODELS 

In  this  section  we  present  the  results  of  a  study  making  use  of  the  solutions 
developed  in  this  article.  The  main  purpose  of  this  study  was  to  evaluate  how 
well  the  classical  square  law  deterministic  and  exponential  interfiring  time 
Lanchester  models  approximate  Models  1.1  and  1.2  (i.e.,  zero-breakpoint-for- 
both-sides  cases)  in  situations  where  we  allow  either  one  or  both  of  the  sides  to 
have  a  gamma(2)  interfiring  time  and  single-shot  kill  probability  which  do  not 
vary  from  round  to  round. 

The  comparisons  were  motivated  by  the  fact  that  it  is  common  in  combat 
models  where  it  is  known  that  the  interfiring  times  are  not  exponential,  to 
assume  they  are  and  use  the  means  of  the  true  distributions.  Thus,  if  p.  is  the 
mean  of  the  true  interfiring  time  distribution,  the  killing  rate  r,  is  taken  to  be 
pip  (done  appropriately  for  each  side)  and  either  the  exponential  Lanchester 
(both  sides  exponentially  distributed)  or  the  deterministic  Lanchester 
differential  equations  are  used  with  the  appropriate  pips  as  the  attrition 
coefficients.  For  a  further  discussion  of  these  matters,  see  references  [2],  [6], 
and  [8]. 

It  should  be  noted  that  in  the  literature  the  exponential  interfiring  time 
model  is  usually  referred  to  as  the  stochastic  Lanchester  model.  In  our  view 
this  is  bad  jargon  because  any  random  interfiring  time  is  stochastic  in  nature. 
We  believe  it  more  accurate  to  modify  the  word  Lanchester  with,  for  example, 
constant,  exponential,  gamma(2),  or  lognormal  in  the  cases  when  the  in- 
lerfiring  times  are  constant,  exponential,  gamma(2),  or  lognormal  random 
variables,  respectively.  Using  the  word  Lanchester,  however,  is  appropriate 
since  the  basic  assumptions  of  independent  firers,  random  selection  of  targets, 
etc.,  is  a  common  thread.  We  also  reserve  the  words  “deterministic  Lanches¬ 
ter”  to  mean  the  classical  differential  equations  of  combat  in  which  there  is  no 
randomness. 

A  perusal  of  the  general  results  of  the  previous  section  and  Appendix  1 
shows  that  the  only  functions  required  to  compute  the  various  state  prob¬ 
abilities  are  the  interkilling  time  density  and  complementary  distribution 
functions  for  sides  A  and  B.  The  generic  versions  of  these  two  functions  for 
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the  constant  round-to-round  exponential  interfiring  time  density  function  and 
single-shot  kill  probability  are  well  known  and  given  by 

=  /SO 

=  0,  otherwise, 
and 

G(/)  =  r>0 

=  0,  otherwise, 

where 

/i  =  mean  interfiring  time, 
p  =  kill  probability. 

For  the  gamma{2)  interfiring  time  random  variable  density  function,  given 
by 

g(/)  =  (4/M-)/e■‘^''‘^  /SO 
=  0,  otherwise, 

we  compute  the  interkilling  time  random  variable  density  function  by  first 
writing  its  distribution  function 

ec 

F(/)=P[T</]=  X  P[r</|N=n]P[N=n] 

n-1 

00 

=  z  HADpq"-', 

J 

where  T  =  interkilling  random  variable,  N  =  round  number  on  which  the  kill 
occurred,  H„{t)  =  n-fold  convolution  of  the  interfiring  time  distribution  H{t) 
with  itself,  p  =  kill  probability,  and  q  =  \  -  p. 

Differentiating  gives  the  density  function 

00 

f(t)=  I,  h„(t)pq''-\  />0 

n  *  1 

=  0,  otherwise. 

Now  the  n-fold  convolution  of  a  gamma(2)  is  a  gamma(2n)  and  is  thus  given 
by 


h„{i)  = 


(2n-l)! 


=  0,  otherwise. 
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so  that  the  above  infinite  series  may  be  written  as 


k  (2n-l)!  ’ 

=  0,  otherwise. 


The  series  on  the  right  is  recognized  as  the  hyperbolic  sine  so  that 
fU)  =  - $\nh[(2lfJL)tq''^l  t>0 

q 

=  0,  otherwise. 

/(f)  may  also  be  written  in  its  exponential  form 

fU)  =  ( >  0 

q 

=  0,  otherwise, 

from  which  we  then  get  by  integration  that 

-(2/m).  ^  ah ^)q"h  -(2/m )q 

«''--27^[Tr77?-TT7^]' 

=  0,  otherwise. 

A  very  important  characteristic  to  note  about  the  two  interkilling  time 
random  variables  is  that  in  the  exponential  interfiring  case  the  interkilling  time 
depends  only  on  the  kill  rate  r-  pifi  but  in  the  gamma(2)  case  both  the  kill 
rate  pip  -  r  and  the  kill  probability  p  (or  mean  interfiring  time  p)  must  be 
specified. 

Our  study  consisted  of  computing  the  relative  difference  that  obtains  when  a 
figure  of  merit  for  a  combat  is  computed  using  the  hypothesized  interkilling 
time  random  variables  for  each  of  the  sides  and  using  either  the  exponential  or 
deterministic  approximation.  Relative  difference  (in  percent)  here  is  defined  as 


where  dc,  -  generic  figure  of  merit  in  the  hypothesized  case,  and  6a  =  the  same 
figure  of  merit  in  either  the  exponential  or  deterministic  approximation  to  the 
hypothesized  case.  Thus,  for  example,  if  both  sides  A  and  B  have  gamma(2) 
interfiring  times  with  parameters  Pa  =  1/H),  Pa  =  1/10  (therefore  Ta  =  1)  and 
Pb  =  2/5,  Ph  -  1/2  (therefore  r;,  =  5/4),  respectively,  then  a  comparison  is 
made  with  sides  A  and  B  both  exponential  (or  deterministic)  with  kill  rates  of 
Ta  =  1  and  rj,  =  5/4,  respectively.  We  also  considered  “mixed”  battles  such  as 
side  A  gamma(2),  side  B  exponential  and  compared  it  with  both  the  exponen- 
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tial-exponential  and  deterministic-deterministic  with  the  appropriate  kill  rates 
for  each  side. 

We  now  discuss  briefly  the  concept  of  parity  and  nonparity  that  we  use  in 
our  subsequent  descriptions.  In  the  deterministic  Lanchester  square  law  model, 
parity  occurs  whenever  rA(ao-  a})  =  raffti)-  b})\  the  battle  goes  to  infinity  and 
neither  side  wins.  We  continue  to  use  the  same  definition  of  parity  in  all 
stochastic  cases  using  r  =  pip..  Strict  parity  is  defined  as  all  parameters  being 
identical  on  both  sides,  i.e.,  the  initial  numbers,  kill  probabilities,  and  in¬ 
terfiring  time  random  variables.  And  when  that  happens  P[A]  =  P[B]  =  1/2.  It 
should  be  noted  that  in  the  exponential-exponential  battles  when  a,)  =  bo  and 
af  =  bf  then  parity  is  equivalent  to  strict  parity. 

The  computations  for  the  state  probabilities  were  made  at  White  Sands 
Missile  Range  on  the  Sperry  1100/82  system.  A  Gaussian  quadrature  tech¬ 
nique  was  used  to  evaluate  all  the  integrals  involved,  not  only  for  the  state 
probabilities  but  the  various  figures  of  merit.  The  particular  quadrature 
technique  we  used  is  described  in  reference  [10]. 

The  accuracy  of  an  integration  over  a  time  interval  [0,  l]  is  a  function  of  the 
number  of  equal  length  segments  the  interval  is  decomposed  into  and  the 
degree  of  the  polynomial  used  for  each  segment.  We  used  three  segments  and 
a  19th-degree  polynomial.  State  probabilities  were  computed  at  30  times 
points  in  the  interval  [0,  /«],  where  ,»  is  essentially  t  =  Thus,  if  one  were  to 
compute  the  state  probabilities  at  u,  the  sum  of  the  transient  probabilities 


Pi2U^)  + P2lit«:)  + Pi2(U)  + PuM  =  0 


and  the  sum  of  the  absorbing  probabilities 


p02(t.)  +  P0l(f-)+Pl0(f«)  +  P2o(t-)  =  1- 


CPU  times  were  the  largest  whenever  gamma(2)-gamma(2)  combats  were 
run,  i.e.,  the  A  side  and  B  side  each  had  gamma(2)  interfiring  times.  And,  for 
these  cases,  the  CPU  times  required  to  compute  all  p,;(t)s,  m^lt),  (r/^(t), 
Weft),  o-aiO  and  the  eight  overall  figures  of  merit  pr^y  ‘T’Td' 
o'[^(°°)].  P[A\  and  P[B]  were  approximately  3  and  2 

hours  for  “reselect  on”  and  “reselect  off,”  respectively. 

Typical  time  varying  characteristics  are  shown  in  Figures  6  and  7.  In  this 
particular  example  we  have  “reselect  on”  (Model  1.1)  and  gamma(2)  in¬ 
terfiring  for  both  sides  A  and  B  with  firing  rates  IIpa  =  10  and  II pa  =  10/9. 
Setting  Pa  =  1/10  and  pa  =  9/10  gives  Ca  =  1  and  rs  =  1.  It  should  be  noted 
that  in  this  case  we  took  i»  =  3.5;  thus  the  sum  of  the  transient  state 
probabilities  at  t«  =  3.5,  in  Figure  6(a),  is  approximately  0  and  the  sum  of  the 
absorbing  state  probabilities  at  =  3.5,  in  Figure  6(b),  is  approximately  1.  The 
30  time  points  at  which  these  state  probabilities  were  computed,  as  shown  in 
Figure  6,  were  selected  so  that  a  Gaussian  quadrature  could  be  used  for 
computing  E[Td]  and  arlTo]  from  Eqs.  (8)  and  (10),  wherein  the  upper  limit 
of  is  replaced  by  /«  =  3.5. 

It  should  be  noted  here  that  the  results  presented  in  Figures  6  and  7  are  for 
a  situation  in  which  parity,  but  not  strict  parity,  obtains;  and  for  the  hypo- 
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TIME 

(b)  Absorbing  Stales 

Figure  6.  State  probabilities  versus  time  for  a  Model  1.1  combat;  A  side  gamma(2) 
with  1/ma=1*K  r^=l;  B  side  gamma(2)  with  IZ/xa  =  I<V9,  Pb  =  9/I0, 

Tr  =  I ;  rescicct  on. 


thesized  interkilling  time  random  variables 

=  P2()(“)  +  Pio('»)  =  P2n(t-)  +  Pio(to.)  =  0.5892, 

P[B]  =  po2(®)  +  Poi{“)  =  P(>2(l-,)  +  pio(f-)  =  0.4106. 

However,  for  the  equivalent  exponential-exponential  combat,  where  r^  =  Za  = 
I,  the  situation  is  one  of  strict  parity  and  P[a]=  P[B]=  1/2.  Figures  8  and  9 
show  the  exponential-exponential  result  corresponding  to  Figures  6  and  7, 
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Figure  7.  Mean  and  standard  deviation  of  number  of  survivors  versus  time  for  a 
Model  1.1  combat;  A  side  gamma(2)  with  =  p.,  =  1/10,  r^  =  l;  B  side 

gamma(2)  with  l/po  =  10/9,  pa  =  9/10,  Tb  =  U  reselect  on. 


respectively.  Note  the  significant  differences  that  occur  in  the  time-varying 
characteristics  between  the  hypothesized  stochastic  and  exponential  models. 

The  differences  between  the  hypothesized  model  and  deterministic 
Lanchester  models  are  even  greater  than  those  obtained  when  the  hypo¬ 
thesized  model  is  compared  to  the  exponential  model.  In  fact,  for  the  case  just 
discussed,  in  the  equivalent  deterministic  case  the  battle  completion  time  is 
infinite  and  neither  side  wins  (since  both  sides  go  to  annihilation).  In  Table  2 
we  present,  for  this  case,  the  values  of  the  eight  overall  battle  figures  of  merit 
for  each  of  the  three  models  and  the  relative  differences  that  are  obtained 
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Figure  8.  State  probabilities  versus  time  for  an  exponential  2-on'2;  A  side  with 
=  1 ;  B  side  with  fg  =  \  . 


when  compared  to  the  hypothesized  model.  Obviously  the  differences  shown 
are  large  and  point  to  the  importance  of  developing  these  models  in  order  to 
obtain  greater  realism. 

It  should  be  noted  here  that  in  our  simulation  studies,  referred  to  in  Section 
1,  we  are  running  much  larger  combats,  for  example,  10()-on-5(),  with  various 
interfiring  time  random  variables  including  the  gamma(2).  We  have  found 
there  that  as  the  battle  size  grows,  i.e.,  as  we  increase  the  initial  numbers 
involved,  the  differences  between  the  hypothesized  model  and  either  the 
exponential  or  deterministic  approximations  gets  larger  percentagewise  also. 
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(b)  Standard  Deviation 

Figure  9.  Mean  and  standard  deviation  of  number  of  survivors  versus  lime  for  an 
exponential  2  of  2;  A  side  with  r^  =  1;  B  side  with  tg  =  1. 


What  is  being  emphasized  here  is  that  we  are  comparing  the  hypothesized  with 
either  the  exponential  or  the  deterministic  approximation  and  not  differences 
between  the  exponential  and  deterministic  models  (although  large  differences 
there  also  exist;  see  reference  [2]). 

In  Appendix  2  is  a  brief  description  of  all  the  cases  run  and  tables  of  relative 
differences  that  were  obtained.  A  perusal  of  those  tables  will  show  relative 
differences,  when  compared  to  the  exponential  model,  as  high  as  44  percent. 
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Table  2.  Overall  figures  of  merit  and  relative  differences  (in  percent)  when  the 
hypothesized  model  is  compared  with  the  exponential  and  deterministic  Lanchester 
models.  For  the  hypothesized  model;  side  A  gamma(2),  l//i^  =  10,  =  I/IO,  =  1; 

side  B  gamma(2),  l/fx^  =  10/9,  pg  =  9/10,  rg  =  1. 


Figure  of 
merit 

Hypothesized 

model 

Exponential 

Deterministic 

Relative 
difference  (%) 
(exponential) 

Relative 
difference  (%) 
(deterministic) 

0.8665 

0.7497 

X 

-1-13.480 

—  X 

0.5359 

0.5568 

0 

-3.900 

100 

£[A(x)] 

0.9919 

0.8337 

0 

-H  5.949 

100 

tr[A(x)] 

0.9016 

0.8974 

0 

+  0.466 

100 

E[B(x)] 

0.6476 

0.8337 

0 

-28.737 

100 

(r(B(3:)] 

0.8376 

0.8974 

0 

-7.139 

100 

P[A] 

0.5892 

0.4997 

0 

4-15.190 

100 

0.4106 

0.4997 

0 

-21.700 

100 

6.  CONCLUSIONS 

Two  versions  of  some  stochastic  homogeneous  two-on-two  combat  models 
have  been  defined  and  state  probabilities  have  been  derived  for  each.  In  the 
first  version,  a  marksman  whose  target  is  killed  resumes  afresh  the  killing 
process  on  a  surviving  target;  in  the  second  version,  the  marksman  whose 
target  is  killed  merely  uses  up  his  remaining  time  to  kill  on  a  surviving  target. 
The  state  probabilities  in  turn  were  used  to  develop  four  time-varying  charac¬ 
teristics  a-Ait),  msit),  crg(t)  and  eight  overall  battle  characteristics 

a-To’  E[A(»)],  or[A(=c)].  £[6(3:)],  a[B(3:)],  P[A].  and  P[B]. 

Comparisons  were  made,  in  terms  of  relative  difference,  with  equivalent 
exponential  and  deterministic  Lanchester  models.  It  was  found  that  both  the 
exponential  and  deterministic  Lanchester  models  are  very  poor  ap¬ 
proximations  of  the  hypothesized  model.  In  fact,  among  the  cases  we  con¬ 
sidered,  there  was  a  relative  difference  of  44  percent  in  the  figure  of  merit 
£[A(3;)]  when  compared  to  the  exponential  model.  And  in  this  study  we 
considered  only  one  interfiring  time  that  was  not  exponential;  namely,  the 
gamma(2).  One  can  reasonably  conjecture  that  larger  relative  differences 
would  surface  when  other  distributions  are  considered.  But  the  point  is  that  we 
have  demonstrated  that  the  exponential  approximation  is  indeed  a  poor  one 
and  that  further  work  must  be  done  to  develop  the  theory  of  small-to- 
moderatc-size  stochastic  combat  models. 


APPENDIX  1 
Model  1.1  (Continued) 

3.  pii(t).  We  now  consider  pnft)  and  write  it  as 

Pi  I ( < )  =  P*iV(  f )  +  P  n’(  t )  +  p  |V(  t )  +  p  n'(  / )■  (13) 


Clearly,  once  we  write  the  first  two  functions  on  the  right-hand  side,  we  may 
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|dvh -  V 


t 

♦i  du  H - u - ^ 

Figure  10.  Definition  of  the  %’ariables  u  and  v  for  the  computation  of  p\"(t) 
and  p,  (i).  X  =  a  B  kills  an  A  in  (t  -  u  -  dv,  t  —  u);  y  =  the  surviving  A  kills  a  B  in 
(t  -  V  —  du,  i  -  u):  u<  V. 


0 


set  down  immediately  the  second  two  by  symmetry.  Consider  now  Figure  10 
above  which  shows  that  a  B  kills  an  A  first  on  the  interval  {i—v  —  dv,  i  —  v), 
the  surviving  A  then  kills  a  B  in  (t—  u  —  du,  t—  u)  with  no  subsequent  killing 
until  after  time  t. 

Now  pn(t)  will  be  written  as 

PiV(t)=  [  dv[  dup\'i{t,  u,  v), 

Jll  Jo 


where  pu\i,  u,  u)  dudv  =  probability  that  both  Bs  are  aiming  at  the  same  A, 
one  of  the  Bs  kills  an  A  in  {t-v-dv,  t-v),  the  surviving  A  kills  a  B  in 
(t-u-du,  t-u),  and  there  is  no  other  killing  until  after  time  t. 
Pn(t,  u,  u)  dudv  may  now  be  written  as  1/2  times  the  product  of 

2/b(»- c)  <Ji>  =  probability  that  one  of  the  two  Bs  kills  an  A  in  the  time  interval 
(t-v-  dv,  I  -  v). 

Ggii  -  v)  =  probability  that  nonkilling  B  has  a  kill  time  >i  -  v. 

Cgiv-  u)  =  probability  that  the  B  which  will  be  killed  by  the  surviving  A  reaims  with 
an  interkill  time  >v-u. 

Gg(v)  -  probability  that  the  B  which  will  survive  reaims  with  an  interkill  time  >  v, 
G^(i-  v)  =  probability  that  killed  A  had  a  time  to  kill  X  -  u, 

/a(<~  u)  du  =  probability  that  the  surviving  A  kills  a  B  in  (( -  u  -  du.  ( -  u), 

and  finally 


Ga(u)  =  the  probability  that  the  killing  A  reaims  with  an  interkill  time  >  u. 
So  we  have 


PtV(f)  =  I  [/eff-  u)CB(f-  t;)GB(t;)GA(t-  u) 

duf^(l- u)Ga{u)Gb(v- u)j.  (14) 

Now  we  consider  pn*(f)  and  write  it  as 
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Given  the  set  of  all  aiming  configurations  in  which  each  B  aims  at  a  different 
/A,  it  is  easy  to  see  that  within  this  set  half  the  time  the  surviving  A  kills  the 
killer  B  and  half  the  time  the  surviving  A  kills  the  nonkiller  B.  So  we  break  up 
Pnft,  u.  u)  into  two  terms.  The  first  term  arises  when  the  surviving  A  kills  the 
killer  and  the  second  term  arises  when  the  surviving  A  kills  the  nonkiller.  Thus 
the  final  result  becomes 


pn’(/.  u,  v)  ciudv  =  ( l/2)"[2/,d  t  -  v)  dv  G„(v  -  u)Gii(l)G,dt  -  v) 

X  u)  du  GaIh)  +  2/3(1-  v)  dv  Gb(v)Gh(i  -  w) 
X  Ga(i-  v)fA(i-  u)  duGA(u)], 

in  which  the  factors  common  to  both  terms  in  the  brackets  are 


l)  dv  =  probability  one  of  the  two  Bs  kills  an  A  in  the  time  interval 
(  /  -  t  -  dv.  I  -  v). 

Oa(i  -  v)  =  probability  that  killed  A  had  a  time  to  kill  >t-v, 
fA(i~  u)  du  =  probability  that  the  surviving  A  kills  a  B  in  (t  -  u  -  du.  1  -  u). 

G^lu)  =  probability  that  the  killing  A  reaims  with  an  interkill  time  >u. 

whereas  the  second  and  third  factors  in  each  term  are  unique  to  the  situation. 
In  the  first  term 

Cb(v  -  u)  =  probability  killer  B  realms  and  has  interkill  time  greater  than  v  -  u  [since  he 
gets  killed  in  the  interval  (t  -  u  -  du.  t  -  u)], 

Gaft)  =  probability  nonkillcr  B  (or  the  surviving  B)  has  an  interkill  time  >t, 

and  in  the  second  term 


Gaft')  =  probability  killer  B  (or  the  surviving  B)  realms  and  has  interkill  time  >v. 
Ga((  -  u)  =  probability  the  nonkiller  B  has  an  interkill  lime  >1-  u  [since  he  gets  killed  in 
the  interval  (1-  u  -  du.  1-  u)]. 


So  finally 

Pn’(f)  =  2  [<Jb(0  I  dv^fB(t-  v)GA(t-  v)  ^  du/A(f- m)Ga(u)Gb(c- u)j 

+  I,  u)GB(u)GA(f-  u)  j  du  GB(f-  u)/A(f-  u)GA(M)j  j. 

(15) 

To  complete  the  computation  for  pii(f)  we  must  write  expressions  for 
Pn(t,  u,  v)  and  p,x(i,  u,  v).  These  two  arise  from  a  transition  from  (2,  1)  to 
(1,  I)  in  which  initially  cither  both  As  are  aiming  at  the  same  B  or  they  are 
not.  In  Figure  11  below  wc  show  that  an  A  kills  a  B  first  in  the  interval 
(I  -  u  -  du,  I  -  u),  the  surviving  B  then  kills  an  A  in  (/  -  u  -  dv,  t  -  v)  with  no 
subsequent  killing  until  after  time  1.  Clearly  we  may  immediately  write  down 
Pi  1(1)  and  pn(f)  by  interchanging  in  equations  (14)  and  (15),  respectively,  A 
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♦I  du  k 


♦I 


■+ 

t 


i  dv  p - V - >\ 


Figure  11.  Definition  of  the  variables  u  and  v  for  the  computations  of  p^ii) 
and  p'n(t)-  y  =  an  A  kills  a  B  in  (i-u-du.  t-u);  j:  =  the  surviving  B  kills  an  A  in 
(t  -  V  -  dv.  V  -  u):  u  >  V. 


with  B  and  u  with  v.  Thus 


Pn’(f)  =  I  du  ^fAii-  u)GA{t-  u)Ga(u)Gb 
x|  f)GB(y)GA(«- t>)|, 


(t-u) 


(16) 


and 


p‘n(f)  =  :^|^GA(f)|  du^fAit-  u)GBit-  u)  ^  du/B(f  -  t'')GB(u)GA{M  -  i;)j 
+  1  du  j^/A(f- m)Ga(m)Gb(J- m)  j  du  Ga(i- u)/B(t- u)GB(i)j  j. 


(17) 


It  should  be  noted  that  in  the  expressions  pVi(f.  u,  v),  /=  1, 2, 3. 4  it  is  not 
always  the  case  that  u  and  v  are  the  backward  recurrence  times,  measured  at 
time  t,  of  the  surviving  A  and  B,  respectively.  Consider  p'n(t.  u,  v)  which 
arises  from  an  initial  configuration  with  each  B  aiming  at  a  different  A.  At 
time  t,  one  of  two  situations  obtain.  Either  the  surviving  B  killed  an  A  or  it 
did  not.  If  it  did  kill  an  A,  its  backward  recurrence  time  is  u;  if  it  did  not  kill 
an  A  its  backward  recurrence  time  is  t. 

Summarizing  to  this  point,  Eqs.  (l)-(3)  and  (13)-(17)  give  all  the  transient 
state  probabilities,  i.e.,  each  of  states  has  in  the  limit  as  f— zero  probability; 
and  we  now  turn  our  attention  to  writing  the  probabilities  associated  with  the 
four  absorbing  states  (0,2),  (2,0),  (0,  1),  and  (1,0). 

d.  P():(f),  p:o(f).  We  now  consider  po2(i)  and  write  it  as 

P()2(f)  =  Poy(t)  +  Po2(l).  (18) 

Consider  Figure  12  below  which  shows  that  a  B  kills  an  A  first  in  the  interval 
(l  -  { -  d{.  i  -  {)  followed  by  another  kill  of  the  surviving  A  by  a  B  in  the 
interval  (t  -  tj  -  drj,  t  -  rj).  Now  po2'(f )  will  be  written  as 

Po2(0=  \  di  {  dTJpI.yfl,  (.  T)), 

Ji)  J{) 
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- - X. - H 

_| - 1  X  I - 1  y  1 - \- 

0  t 

- HdTih - ^ - H 

Figure  ^12.  Definition  of  the  variables  f  and  rj  for  the  computations  of  pl^^i) 
and  Posit)-  Jt=a  B  kills  an  A  in  i-f);  y  =  a  B  kills  another  A  in 

(/-  77  -  d-n.  t-T}):  {>r). 


where  po'sV,  17)  dCdr)  =  probability  that  both  Bs  are  aiming  at  the  same  A, 
one  of  the  Bs  kills  an  A  in  ^  -  t  -  and  then  one  of  the  Bs  kills  the 
surviving  A  in  (f-  tj-  dr),  t-  77).  77)  d^dr)  may  now  be  written  as  1/2 

times  the  product  of 

2/b(<“  f)  di  =  probability  that  one  of  the  two  Bs  kills  an  A  in  the  time  interval 

GbU  ~  0  =  probability  that  nonkilling  B  has  a  kill  time  >1-  C, 

-  faU  “  7))  47J  =  probability,  after  both  Bs  reaim,  that  a  B  kills  the  surviving  A  in  the  time 
interval  (i-  rj-  dr),  t)), 

Gb((~  tj)  =  probability  that  nonkilling  B  has  an  interkill  time  >{  -  r), 

~  probability  that  the  first  A  that  was  killed  had  a  kill  time  >1-  {, 

and  finally  v)  -  probability  that  the  second  A  that  was  killed  had  a  kill 

time  >1-7).  So  we  have 

p;,y(/) = 2  di:  [/s(t  -  ooso  -  ogaU-  0 

I  dvfsiC-  v)Gb(C-  v)GA{t-  v)^.  (19) 

Now  consider  pozit)  and  write  it  as 

pilaff)  =  [  dC  \  dr)pokt,Lv)- 

Jo  Jo 

Here  we  must  break  up  pJzil,  rj)  into  two  terms.  The  first  term  arises  when 
the  same  B  kills  both  As  and  the  second  term  arises  when  each  B  kills  an  A. 
The  final  result  becomes 

Pozft,  C,  ■h)  d(dr)  =  \l2[2foU- 0  dCfeiC-v)  dr]  GbO  -  r))GA{t  ~  OGa^i  -  17) 
+  2f„(t-  0  dif„(t-  77)  dr]  Gb((-  r])GA((-  ^)Ga(<-  t?)] 

in  which  the  factors  common  to  both  terms  in  the  brackets  arc. 


2 fn( (“{)  “  probability  one  of  the  two  Bs  kills  the  first  A  in  the  time  interval 

U-i-d{.  »-{), 
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GaIi-  0  =  probability  that  the  first  killed  A  had  a  lime  to  kill  >  /  -  f, 

Ga(i  ~  v)-  probability  that  the  second  killed  A  had  a  time  to  kill  >  t  -  tj, 

whereas  the  factors  not  common  are  unique  to  the  situation.  In  the  first  term 

r|)  dri  =  probability  that  the  killing  8  reaims  and  kills  the  surviving  A  in  the  lime 
interval  (l-  rj-  drj.  t  -  77), 

CbU~  n)  =  probability  that  the  B  who  kills  no  As  has  a  kill  time  X  -  tj, 
and  in  the  second  term 

/a(t  -  ij)  dtj  =  probability  that  the  other  B  (who  did  not  have  to  reaim)  kills  the  surviving 
A  in  the  interval  (/  -  tj  -  dr),  t  -  tj), 

GbH-  tj)  =  probability  that  the  first  B  who  killed  an  A  (and  had  to  reaim)  has  an 
interkill  time  >  f  -  t). 

So,  finally, 

P02(0  =  £  d(l^/s(t-0GA(t-0  1^'  dvGs(l-v)fBU-v)GAU-v)] 

+  £  j'  dvGsU-  T))/e(f-  T,)G^(t~  T,)]. 

(20) 

Similarly,  write 

P2o{t)  =  P2o(r)  +  P2n(t).  (21) 

If  we  suppose  that  the  first  B  is  killed  in  (f  -  ^  -  t-  0  and  the  second  B  is 
killed  in  (t-  r}-  dr),  t-  rj),  then  by  simply  interchanging  A  with  B  in  Eqs. 
(19)  and  (20)  we  may  get  p*2o(i)  and  p^zoit).  Thus 

pS(«)  =  2^^dc[fAt-  OGaH-  OGsd-O 

x£  dv/Ad  -  v)GaU-  v)GB(t~  v)]  '22) 

and 

dUfAit-OGBit-O  dr,  Ga(/-  t,)/a(C-  T,)GB(t-  T,)] 

+  £  d(^Ut  -  ^)G„((  -  0  j  ^  dr,  GaU-  v)fA(l-  v)Gb(  t  -  t,)]  ■ 

(23) 

5,  P()i(t),  Pto(<).  We  next  consider  poift)  and  write  it  as 

P.||(<)  =  Pui*(l)  +  p!n(f)  +  p,'h  (f)  +  Pin(t). 


(24) 
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To  develop  expressions  for  the  terms  on  the  right-hand  side  of  Eq.  (24)  it  is 
very  convenient  to  work  with  the  backward  recurrence  time  and  its  associated 
instantaneous  kill  rate.  For  the  first  term,  recall  both  Bs  are  initially  aiming  at 
the  same  A  and  refer  back  to  Figure  10  and  consider  the  following  sequence 
of  events: 


1.  One  of  the  Bs  killed  an  A  in  the  time  interval  (f  -  u  -  dv,  t  ~  v). 

2.  Then  both  Bs  reaimed  ai.d  one  of  the  Bs  was  subsequently  killed  by  the  surviving  A 

in  the  time  interval  (/  -  n  -  du.  t-  u). 

Therefore,  at  time  t  the  surviving  A  has  a  backward  recurrence  time  of  u  with 
instantaneous  kill  rate  TaIu)  and  the  surviving  B  has  a  backward  recurrence 
time  of  V  with  instantaneous  kill  rate  raiv). 

Consider  now  the  time  interval  (t.  /  -(-A)  and,  in  the  usual  fashion,  write  an 
expression  for  pI.|’(/  +  A)  retaining  only  the  first  order  terms  in  A.  Thus  the 
probability  of  being  in  state  (0,  1)"’  at  time  f  +  A  is  equal  to  the  sum  of  two 
probabilities,  namely,  (1)  the  probability  of  being  in  state  (0,  l)‘'*at  time  t 
times  the  probability  of  remaining  there  in  (f,  t-t- A)  [which  is  one  since  (0,  1)“’ 
is  an  absorbing  state],  and  (2)  the  probability  u.  v)  du  dv.  of  being  in  the 
state  {1,  1)“’  at  time  t  with  backward  recurrence  times  of  u  and  v  for  survivors 
on  the  A  side  and  B  side,  respectively,  times  the  probability,  rB(u)A(l  - 
rA{u)A),  that  the  B-side  survivor  kills  the  A  and  A  fails  to  kill.  Thus,  we  get 
taking  into  account  all  (u,  v)  pairs. 


pS)V(t  +  A)  =  poV{t)+  [  dv  [  dupu(t,  M,  v)rB{v)A{]  -  rA(u)A). 

Jo  Jo 


Rearranging  terms,  dividing  by  A,  and  letting  A-»0.  we  have  that 


dp\u(i)/di=\  dv  \  du pu{t.  u,  v)rBiv). 
Ji)  Jo 

Using  the  initial  condition  poV(O)  =  0,  we  may  write 

PoV(f)=  [  dC  {  dv  {  dup\\\c.  u,v)rB{.v). 

Jo  Jo  Jo 


Finally,  after  substituting  in  for  rofu)  =  faMIGsiv)  and  u,  u)  [see  Eq. 

(14)],  we  get 


PoV(f)=  £  d(  I'  dv  ^fBiv)fB((~  v)GbU-  v)GM-v) 

x|  duGidv  -  u)fA{i  -  u)Ga(u)^.  (25) 

Proceeding  in  the  above  manner,  i.e.,  making  use  of  the  backward  recur¬ 
rence  times  and  the  instantaneous  kill  rate  associated  with  each,  it  is  easy  to 
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get  the  next  three  equations,  namely, 


Po.V)  =  ^  {  dC  [/b(^) dv  [/b(^  -  v)GAC  -  v) 
x|  du  GbIu- u)/a(^- u)GA(M)j  j 
+  di  dv^feM/sU -  L')GAi-  v) 

dM  Gb(^- u)/a(^- w)GA(u)j|, 

Pm(0  =  |  dC^^  tin  |^/a(^- w)GA(f- m)Ga{w)Gb(^- «) 
x|  du/BC^- n)/B(u)GA(u- u)j. 


(26) 


(27) 


and 


pin  (0  =  5  {|/f  [  Ga(C)  du  [UC-  u)  GsU  -  u) 
xj  dt)/s(u)/B(^- 1))Ga(«- I>)jj 
+  j^u  [/a(^-  «)Ga(m)Gb(^-  n) 

X  dv  GM-  i^)/b(^-  i))/B(t')]).  (28) 


We  next  write  the  right-hand  terms  of 

Pio(0  =  p"o(0  +  p'io(f)  +  p‘io(t)  +  p'io(0-  (29) 

Each  of  these  are  easily  written  by  taking  the  appropriate  pln(i),  i  =  1,  2,  3,  4 
and  interchanging  A  with  B  and  u  with  u.  Table  3  below  lists  these 
identifications.  Thus  to  get  pio(f)  interchange  A  and  B  and  u  with  v  in  poift), 
etc.  These  identifications  result  in  the  formulas 

p'.«(/)  =  £  di  £  dv  [fBU-v)GB(C- v)Gb(v)GM{- v) 

x('’dufM-u)fA{u)Gs(v-u}],  (30) 
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p'.o(f)  =  ^{£  dv^fB{C-v)GAC-v) 

x|^  duUu)fAC-u)GB(v-u)J^ 

X  t/uGB(f-w)/^(f-u)/A(u)]},  (31) 

p'o(t)  =  dC  £  dll  [/a{m)/a(^-  u.)Ga{C-u)GbU-u) 

i>)/b(^- u)GB(t;)j,  (32) 

and 


Pio(0  =  \  {£  dC  |^/^(^)  du  [/a{^-  u)Gb{C-u) 

ft  U 

xj^  dvGAu~v)fB{C-v)GB{v)J^ 

A  ^ 

x|^  t/t  G^(^-t;)/B(f-t;)GB(t;)]}.  (33) 


Table  3.  Identifications  used  to  get  p,o(0  from  Podt). 


(1. 1)^(0, 1) 

(1. 1)^(1, 0) 

Pol(t) 

Initial  aiming 
configuration 
and  winning  side 

(i),  V 

PwU) 

Initial  aiming 
configuration 
and  winning  side 

PoV(t) 

Both  Bs  aiming 
at  same  A  and  B 
side  wins 

p',»(t) 

Both  As  aiming 
at  same  B  and  A 
side  wins 

Po’ft) 

Each  B  aiming  at 
a  different  A  and 

B  side  wins 

p',o{t) 

Each  A  aiming  at 
a  different  B  and 

A  side  wins 

P«,V) 

Both  As  aiming 
at  same  B  and  B 
side  wins 

Pio(f) 

Both  Bs  aiming  at 
the  same  A  and  A 
side  wins 

Poi’(t) 

Each  A  aiming  at 
a  different  B  and 

B  side  wins 

(2).  . 
PwU) 

Each  B  aiming  at 
a  different  A  and 

A  side  wins 
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Model  1.2  (ao  =  2.  af  =  0;  bo  =  2,  b,  =  0;  “Reselect  Off”) 

The  reader  is  reminded  that  in  Version  2  of  Model  1  (Model  1.2),  i.e., 
“resclect  off,”  if  a  marksman’s  target  is  killed  by  the  other  member  on  the  side, 
his  remaining  time  to  a  firing  (or  a  killing)  is  carried  over  to  the  survivor  (if 
there  is  one).  To  develop  the  state  probabilities  for  this  model  it  is  no  longer 
necessary  to  break  down  the  analysis  into  terms  of  the  initial  aiming 
configurations,  i.e.,  both  Bs  aiming  at  the  same  A,  etc.  However,  we  still  do 
decompose  state  (1,  1)  into  states  (1,  1)*”,  (1.  1)*‘\  (I,  1)‘^*,  and  (1,  I)'"*'  and 
these  are  defined  in  Figure  13  below. 

/.  p22i0-  As  in  Model  1.1,  we  have 

P22(n  =  (G^(t))-(GB(t))-.  (34) 

2.  pi2(0,  Piiit)-  Referring  back  to  Figure  5  we  define  pizit,  17)^7?  =  one 
of  the  Bs  kills  an  A  in  (t  -  p  -  dp.  t  -  p),  and  there  are  no  killings  until 
beyond  r.  Now 

Pi2(t.  v)  dn  =  -/e(t-  9)  dp  Gs(t?)Ga(/-  p)GAit), 


where 

2/B(f~  7)  4»7  =  probability  that  one  of  the  two  Bs  kills  an  A  in  the  time  interval 

(t  -  t)  -  di).  i  -  Tj). 

Gb(7)  =  probability  that  the  B  that  killed  reaimed  and  has  an  interkill  time  >t), 

GbU)  =  probability  that  the  nonkiller  B,  since  he  just  carries  over  his  remaining 
time  to  kill  if  he  is  aiming  at  the  killed  A,  has  a  time  to  kill  >1, 

Ca(‘-  7)  =  probability  that  killed  A  had  a  time  to  kill  >1-  p, 

and  finally  G^(f)  =  probability  that  surviving  A  had  a  time  to  kill  >  t.  Thus  we 
get 


p,2(t)=[  dppnil.  p)  =  2Gb(i)Ga(i)  I  dpfsic- p)Gb(p)Ga(i- p). 

Jo  Jo 


(35) 


(1.2) 


(1.1)<') 


(1,1)(2) 


surviving  B  in  stale  (1,1)  killed  A 


surviving  B  in  state  (1,1)  did  not  kill  A 


(2.1) 


Figure  13. 

off”  model. 


(1.1) 


(3) 


(1.1)H) 

Definition  of  the  states  1. 


surviving  A  in  stale  (1,1)  killed  B 

surviving  A  in  stale  (1,1)  did  not  kill  B 

))"’,  (1,  1)"’,  (1,  1)'”,  and  (1.  1)'"’  in  the  “rcselcct 
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Again,  by  an  interchange  of  A  with  B,  we  get 

P2,(/)  =  2G^(t)GB(0  f  diUil- OGM) Gait- 0-  (36) 

Jo 

3.  pii(r).  As  described  above  we  shall  write  Pi\{t)  as 

Pi  i(»)  =  Pn’(f)  +  Pn’(t)  +  Pn’(0  +  PuU)-  (37) 

For  the  p*iV(0  and  pu*(<)  computations  refer  to  Figure  10.  After  the  first  kill  of 
an  A  by  a  B,  it  is  clear  from  the  independence  assumption  and  homogeneity  of 
the  marksman  on  each  side  that  with  equal  likelihood  the  surviving  A  is  either 
aiming  at  the  B  that  killed  or  the  B  that  did  not  kill.  Thus  define 
Pi  i(/.  w,  v)  dudv  =  probability  that  one  of  the  Bs  killed  an  A  in  (r-  t;  -  dv,  t- 
u),  the  surviving  A  kills  the  nonkilling  J5  in  (/  -  u  -  du,  t-  u),  and  there  is  no 
subsequent  killing  until  after  time  /;  u,  v)  dudv  may  now  be  written  as 

1/2  times  the  product  of 


IfsU-  u)  dv  =  probability  that  one  of  the  two  Bs  kills  an  A  in  the  time  interval  (f  -  v- 
dv,  t-  v). 

Ggiv)  =  probability  that  the  killing  B  reaims  and  has  an  interkill  time  >  v  (since  this 
is  the  B  that  will  survive  on  the  B  side). 

Gfifi-  u)  =  probability  that  the  nonkilling  B,  which  will  be  killed  by  surviving  A,  has  a 
time  to  kill  >t  -  u  (this  is  true  whether  he  were  initially  aiming  at  the  A 
that  was  killed  or  not,  since  if  he  were  aiming  at  the  killed  A,  he  merely 
transfers  the  remaining  time  to  kill  to  the  surviving  A), 

GaU-  ti)  =  probability  that  killed  A  had  a  time  to  kill  >t  -  ti, 
fA(i~  M)  du  =  probability  that  the  surviving  A  kills  a  B  in  (i-  u-du,  t-  u), 

and  finally  Ga(m)  =  probability  that  the  killing  A  reaims  with  an  interkill  time 
>«.  So  we  have 

pi‘i’(f)  =  |  duj^/eff- u)GB(i>)GA(f- u)  duGait-  u)fA{t-  u)Ga{u)^- 

(38) 


In  a  similar  fashion  we  define  p*n  (f,  u,  v)  du  dv  =  probability  that  one  of  the  Bs 
killed  an  A  in  (i  -  v  -  dv,  t  -  v),  the  surviving  A  kills  the  killing  B  in 
(i-  u-  du,  t  -  u),  and  there  is  no  subsequent  killing  until  after  time  t,  and  get 

Pn(f,  u,  v)  dudv=  fait  -  v)  dv  Gaiv  -  u)GBit)GAit-  v)fA{t-  u)  du  Ga(u), 
with 


pmV) 


=  Gait)  dv^fait-  v)GAit-  v)  ^  d«  Gb(i;  -  u)/a(<  -  «)Ga(u) j  • 


(39) 


To  write  expressions  for  p^n(f)  and  pn'(f)  see  Figure  11  to 


recall  the 
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definitions  of  u  and  v  in  the  transition  (2, 1 )-+(!,  1).  As  before,  we  can,  by 
interchanging  A  with  B,  u  with  i;  in  Eqs.  (38)  and  (39)  set  down 

Pn'(f)  =  |  du^fAit-  u)GAiu)GB{t-  u)  j  dt)  GA(t- u)/B(t- t;)GA(u)j, 

(40) 

and 


Pn(t)  =  GA(t)  I  du  ^/aO-  u)GB{t-  u)  j  dv  Ga(u  -  v)fBit-  i))GB(u)j, 

(41) 

respectively. 

4.  poiU),  PioiO-  We  now  write 

Polit)  =  P02(t)  +  P02(»),  (42) 

where  poVft)  is  the  term  that  arises  when  the  same  B  kills  both  As  and  poW) 
arises  when  each  B  kills  an  A.  Referring  to  Figure  12  for  the  definitions  of  { 
and  7]  we  may  write 

Po2(t)  =  2  d{l^fB(t-OGA(t-C)  dr,fBU-v)GB(t-v)GA(t-v)], 

(43) 

and 


p^02\ 0  =  2  ly^l^Ut-OGAd-oj^  dvGsU-  v)/b(i  -  v) Ga( I  -  t,)]  . 

(44) 

Similarly,  write 

P2o(t)  =  P2o(t)  +  P2o(0-  (45) 

If  we  suppose  the  first  B  is  killed  in  (l  -  d^,  t  -  and  the  second  B  in 

(t  -  T)  -  drf,  dj}),  then  by  simply  interchanging  A  with  B  in  Eqs.  (34)  and  (35) 
we  get 

P2o(t)  =  2  d^j^Ut-OGBU-O  dr,fAU-v)GA(t-v)GB(l-  t,)], 

(46) 

and 


p‘2o(0  =  2  dC  [fA{t-0  GbO-O  dTj  GaU- v)fA(t-v)GB(l-v)l 


(47) 


respectively. 
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5.  Pw(t)-  We  write 

Poi(0  =  Poi’(»)  +  po!(t)  +  Po,’(0  +  poi’(r)  (48) 

and  use.  as  in  Model  1.1  instantaneous  kill  rates,  each  of  which  depends  on  the 
appropriate  backward  recurrence  time,  to  develop  the  expressions  for  the 
terms  on  the  right-hand  side  of  Eq.  (48).  Thus  in  the  case  of  p^'/ff),  which 
arises  from  the  sequence  (2, 2)-^  (1, 2)->  (1,  l)->  (0,  1),  we  take  pu(t,  u,  u),  in 
which  u  and  v  are  the  backward  recurrence  times,  measured  at  time  t,  of  the 
surviving  A  and  B,  respectively,  and  in  the  usual  fashion  write 

poV(f  +  A)  =  poV(0+  [  dv  \  dupuit,  u,  v)rB(v)A(l  -  u(u)S). 

Jo  Jo 

After  rearranging  terms,  dividing  by  A,  letting  A-»0,  and  using  the  initial 
condition  poi(O)  =  0,  we  may  write 

poV(r)  =  dv  [fBU-v)fB(v)GAi-v) 

x|  du  Gb(^- u)/a(^- M)GA(u)jj.  (49) 

In  a  similar  fashion  we  also  get 

Po.’(t)  =  dc  dv  [/B(f-  i^)Ga(^-  u) 

x£  dM/A(^-u)GA(u)GB(i;-u)]],  (50) 

pt'f(0  =  [  d^[|'  du^fAC-u)GAu)GBU-u) 

x£  dvGM-v)fB(C-v)fB{v)^],  (51) 

and 

P<n(t)  =  |'  dc\^GA0  du]^U(-u)GBU-u) 

x|  ii)/s(i>)Ga(u  -  u)jj.  (52) 

Again,  as  in  Model  1.1,  we  may  interchange  A  and  B  and  u  with  v  in  the 
appropriate  pln(i),  i  -  1,2,  3, 4  to  get  each  of  right-hand  members  in 

Pio(t)  =  Pl()(t)  +  Plo(0  +  PlC)(t)  +  Plo(0 


(53) 
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plo(t)  =  £  [  j'  dv  [fB((-v)GB{v)GM- v) 

x|  d£/GB(^- «)/A(w)jj,  (54) 

p'.o(0  =  [Gb(^)  dv  [/b(^-  v)GA(-  v) 

^/«/a(^-m)/a(u)Gb(u- u)jj,  (55) 

p‘.o(0  =  dC  [  du  [/a(^-  u)/a(m)Gb(C-  «) 

x|  duGA(C- u)/b(^- t))GB(tJ)jj,  (56) 

p'.oV)  =  Ij  di  [/a(^)  du  [fM-u)GB{C-u) 

i^)Gb(u)Ga(u- D)jj.  (57) 


Model  2.1  (So  =  2,  a,  =  1;  bo  =  2,  6i  =  0;  “Reselect  On”) 

The  transient  states  are  given  by 

Piiil)  ~  (Ga(0)^(Gb(/))^, 

P2l(f)=GB(f)  f  dT)/A(f-  T?)GA(f-  Tj)(GA(T)))^GB(t-  Tj) 
Jo 

+  Gb(/)Ga(/)  {'  dvfAU-v)GA{v)GB{t-v)- 
Ji) 

The  absorbing  states  are  given  by 

P,2(»)  =  2  {'  diiGAiOfCBiOfaU), 

Jo 

Pn(0=f  drj[  dM/A(T7- u)Ga(t|- u)(Ga(u))^Gb(t}- u)/b(t}) 
Jo  Jo 

+  \  drj  {  d«/A(T7- «)Gs(t7- u)Ga(m)Ga(i7)/b(7)), 

Jo  Jo 


(58) 


(59) 


(60) 


(61) 
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p:o(0  =  2  £  -  OGbU  -  0  driUi- v) 

X  Ga(C-  v)(jb(i-  T?)j 

+  1  dr]GA(l-v) 

xfAU-v)GBO~v)j 

+  £  d^^fAd-  OGbU-C)  df)  GaU-v) 

^/a(<~  17)Gb(*- i7)j.  (62) 

Model  2.2  (ao  =  2,  a/  =  1;  bo  =  2,  6#  =  0;  “Reselect  Off’) 

The  transient  states  are  given  by 

P22d)  =  (GA{t))-{GB[t))\  (63) 

P2i(f)  =  2GA(t)GB(f)  [  dCfA(t-0GA(0GB(l-0-  (64) 

Jo 

The  absorbing  stales  are  given  by 

p,2(0  =  2f  (GA(0fGB(0/B(()d(.  (65) 

Jo 

pi,(0  =  2[  drf  {  du/AlT)  -  u)Ga(u)Ga(t})Gb(tj- «)/b(i7).  (66) 

Jo  Jo 

P2()(/)  =  2|  d^^fAd-OGad-O^  drifAU- v)GAd- v)GBd- v)^ 

+  2|^  dC^fAd-OGad-C) 

x|  dTJ  Ga(^- 1?)/A(f-  17)GB(t-  T})j.  (67) 

Model  3  (ao  =  2,  ar  =  1;  Ih>  =  2,  Ar  =  1;  “Reselect"  not  Material) 

P22(f)  =  (GA(^))'(GB(/))^  (68) 

while  the  absorbing  state  probabilities  are 

p,2(t)  =  2f  {GA(C)?GBU)fB(C)d(,  (69) 

Jo 
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P2,(t)  =  2[  (GB(0}-GM)fAU)d{.  (70) 

Jo 


APPENDIX  2 

As  discussed  in  Section  4  we  define  parity  to  mean  =  rabo  (since  we  are 
assuming  here  that  af  =  bf  =  0)  so  that  if  ao  =  bo  parity  is  equivalent  to  =  Tb 
and  nonparity  is  equivalent  to  Ta  rs.  Table  4  shows  all  the  gamma(2)- 
gamma(2)  combats  that  we  considered  in  the  nonparity  situation.  Thus  we 
have  here  nine  cases  for  “reselect  on”  and  nine  cases  for  “reselect  off”  for  a 
total  of  18  different  cases. 

We  also  ran  the  same  probabilities  displayed  in  Table  4  for  the  parity  case 
ta  =  rfl  =  1.  Again  there  were  a  total  18  different  cases  for  the  hypothesized 
model  (including  “reselect  on”  and  “reselect  off").  It  should  be  noted  that  in 
the  hypothesized  case  strict  parity  obtains  only  when  Pa  =  Pb!  also,  the  results 
for,  say,  pa  =  9/10  and  Pb=1/10  can  be  obtained  by  interchanging  the 
outcomes  from  the  Pa  =  1/10  and  pe  =  9/10  battle.  The  “mixed”  battles,  one 
side  gamma(2)  the  other  side  exponential,  included: 

1  A  side  gamma(2).  Ca'-I-  />a  ='/!().  1/2.  9/10;  B  side  expDnential.  rn=1.25; 
nunparity.  Run  with  •Teselecl  on"  and  •reseleei  off  for  side  A. 

2.  A  side  exponential.  rA=\\  H  side  gammal2).  r„=1.25,  p»=l/IO.  1/2.  9/10; 
nonparity.  Run  with  "reselect  on"  and  "reselect  off"  for  side  B. 

.t,  A  side  i:amma(2).  r^  =  1.  Pa  =  I/O*-  1/2.  9/10;  B  side  exponential.  =  I;  parity. 

Run  with  "reselect  on"  and  "reselect  off"  for  side  A. 

4.  A  side  exponential,  r a  =  I:  B  side  gamma(2).  Tb  =  I.  pa  =  1/10.  1/2.  9/10;  parity. 

Run  with  "reselect  on"  and  "reselect  i>ff"  for  side  B.  iNitte  that  the  results  here  would 
he  obtained  by  just  interchanging  the  A  side  and  B  side  results  in  .3,1 

To  make  the  relative  difference  calculations  discussed  earlier,  only  two 
exponential-exponential  and  two  deterministic-deterministic  battles  need  be 
considered;  namely,  rA=1.00,  re  =125  and  rA=1.00,  rB  =  1.00  for  each 
pair. 

Tables  5-14  which  follow,  present  all  relative  differences  obtained  when  the 
hypothesized  model  is  compared  with  the  exponential  model.  We  have  also 
included,  for  the  sake  of  completeness,  some  relative  difference  results  when 
the  hypothesized  model  is  compared  with  the  deterministic  model.  These  are 
shown  in  Tables  15  and  16.  Also  in  Tables  7,  8,  13,  and  14,  where  there  are 
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results  shown  for  the  exponential-exponential  case  with  parity,  i.e.,  =  /"b  = 

1,  the  win  probabilities  are  entered  as  P[ A]  =  F[B]  =  0.4997.  This  is  the 
answer  obtained  using  the  Gaussian  quadrature  (we  know  that  we  must  have 
P[A]=  P[B]=  1/2). 
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